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We show that isotropic component of circular polarization of stochastic gravitational wave back- 
ground can be explored by breaking two dimensional configuration of multiple laser interferometers 
for correlation analysis. By appropriately selecting orbital parameters for the proposed BBO mis- 
sion, the circular polarization degree II can be measured down to II ~ 0.08(fi G w /10" 15 )~ 1 (SiV J R/5) 
with slightly (~ 10%) sacrificing the detection limit for the total intensity Qgw compared to the 
standard plane symmetric configuration. This might allow us to detect signature of parity violation 
in the very early universe. 

INTRODUCTION 

<N ■ 

Due to extreme penetrating power of gravitational waves, observation of the waves may enable us to study the very 
, early universe in a way that cannot be achieved with other methods Q. Recently, follow-on missions to the Laser 
Interferometer Space Antenna (LISA) Q have been actively discussed to directly detect stochastic gravitational wave 
, background from the early universe around ~ 1Hz. The primary target for the proposed missions, such as, the Big 
Bang Observer (BBO) 0\ and DECIGO 0], is the gravitational wave background produced at inflation. With various 
observational supports for inflation, existence of the background is plausible, while it is currently difficult to predict 
' its amplitude. 

Meanwhile, considering the facts that observation of gravitational waves will be a truly new frontier of cosmology 
' and our understanding of physics is limited at very high energy scale, it would be quite meaningful to prepare flexibly to 
various observational possibilities. Actually, detection of something unexpected with odd properties will be generally 
^\ ' more exciting than confirmation of something widely anticipated. For this purpose it is preferable that we can study 
gravitational wave background beyond its simple spectral information, and model independent approach would be 
effective with regard to the sources of the background. 

One of fundamental as well as interesting properties of the background is its circular polarization. Circular po- 
' larization characterizes asymmetry of amplitudes of right- and left-handed waves, and its generation is expected to 
be closely related to parity violation (see e.g. 0). In a recent paper 6] it was shown that LISA can measure the 
dipolc (I = 1) and octupole (I — 3) anisotropic patterns of circular polarization of the background in a relatively clean 
, manner, but, at the same time, LISA cannot capture its monopole (/ = 0) mode due to a symmetric reason. Since 
our observed universe is highly homogeneous and isotropic at cosmological scales, it would be essential to have sen- 
sitivity to the monopole mode of circular polarization of cosmological background. The proposed missions like BBO 
or DECIGO are planed to use multiple sets of interferometers to perform correlation analysis to measure the total 
intensity flaw of cosmological background by beating out detector noises with a long-term signal integration. The 
standard configuration of these missions is to put the multiple interferometers on a same plane. This is advantageous 
to get a good sensitivity to the total intensity Qgw, as a larger spatial separation between interferometers results 
in reducing their overlapped responses to gravitational waves [jj. However, with this plane symmetric configuration, 
we are totally blind to the monopole mode of circular polarization, as in the case of LISA. This means that even if 
the isotropic background is circularly polarized by 100%, we will not be able to discriminate its extreme nature. In 
this paper we study how well we can detect the monopole mode of circular polarization by breaking symmetry of the 
detector configuration. 
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CIRCULAR POLARIZATION 



In the transverse-traceless gauge, the metric perturbation due to gravitational waves is expressed by superposition 
of plane waves as follows: 

h ab (t,x)= J2 I" df I dnhp^ny^-n-Ve^n), (f) 

P=+ x 

where S 2 is the unit sphere for the angular integral, the unite vector n = (sin 9 cos (f>, sin 9 sin <fi, cos 9) is propagation 
direction of the waves, and e+ fc and e* b are the bases for transverse-traceless tensors. We fix them as e + = eg eg) eg — 
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® and e x = eg ® + ® eg where two unit vectors eg and are denned in a fixed spherical coordinate 
system, as usual. The frequency dependence is easily resolved by Fourier transformation, and we omit its explicit 
dependence for notational simplicity, unless we need to keep it. We decompose the covariance matrix (hpi(n)h P2 (n')) 
(PI, P2 = +, x) for two polarization modes as 



6 drc (n-n r ) f I( n ) + Q(n) U(n) - iV(n) 
4tt \U(n)+iV(n) I(n) - Q(n) 



(2) 



where I,Q,U and V are the Stokes parameters and are real by definition. The parameters Q and U are related to 
linear polarization, and their combinations Q ± iU have spin ±4 and are expanded in terms of the spin-weighted 
spherical harmonics ± 4 Y; m (n) defined for I > 4 [fj. Therefore, they do not have modes naturally corresponding to 
the monopole pattern. The parameter / represents the total intensity, while the parameter V shows asymmetry of 
amplitudes of right- and left-handed waves. These meanings become apparent with using the circular polarization 
bases br^l = (e+ ± £e x )/v2 for which the expansion coefficients become hn t L = (h+ =F i/i x )/v2- The parameters 
/ and V can be expressed as I(n) = (\h R \ 2 + \h L \ 2 ) /2 = (\h+\ 2 + \h x \ 2 )/2 and V(n) = (\h R \ 2 - \h L \ 2 ) /2 = 
i(h + h* x — h±h x ) /2. They are spin quantities and have monopole modes 7 o and Voo for spherical harmonic 
expansions I(n) — Ylim IlrnYi m (n) and V(n) = J2im VlmXlm( n )- 111 this paper we mainly discuss how to capture the 
monopole mode of circular polarization signal Voo with laser interferometers in space. 

To begin with, we summarize basic aspects of LISA [j|. LISA is formed by three spacecrafts that nearly keep a 
regular triangle with its side length L = 5 x 10 6 km. The geometric barycenter of the triangle moves on a circular 
orbit around the Sun with radius ~ 1AU. The detector plane made by the triangle inclines to the orbital plane of 
the barycenter by ~ 60° with changing its orientation. The triangle also rotates annually on the detector plane (the 
so-called "cartwheel motions"). From its six one-way relative frequency fluctuations of the laser light, we can make 
three Time-Delay-Interferometer (TDI) variables {A, E, T} that cancel the laser frequency noises The detector 
noises between these three variables are not correlated due to a symmetry of the data streams, and two modes {A, E} 
have similar noise spectra 0- At low frequency regime / = 2irf /L <C 1, the responses of two modes {A,E} to 
gravitational waves can be regarded as those of two simple L-shaped interferometers that measure differences of two 
arm- lengths caused by passing gravitational waves, as shown in figure 1. The T-mode is less sensitive to gravitational 
waves at / < 1 and is not important in this paper (see e.g. 0,^3)- Since the monopole modes loo an d Voo are our 
main concern and they are invariant under rotation of a coordinate system, we only use the coordinate system XYZ 
fixed to the single system Bl as shown in figure 1. These basic aspects for LISA are essentially same for BBO 
But BBO is planed to have a smaller arm-length L = 5 x 10 4 km (/ = 1 corresponding to 0.95Hz) and use multiple 
systems (triangles) not only the first system Bl [j| (see also 0,0|). The responses of A and E modes to gravitational 
waves are written as 

{A,E}= f dn Y, h P (n)F( AE} (nJ), (3) 

52 P=+,x 

and the pattern functions F^(n,f) for A mode are expressed as F^(n 1 f) = |(1 + cos 2 9) sin(2</>) + 0(f) and 
F^(n, f) = cos 9 cos(2</>) + 0(f) (see e.g. 0,0]). The functions F^(n, /) for E mode is given by replacing with 
4> + 7r/4. Here we multiplied an appropriate common factor proportional to some powers of / so that the pattern 
functions become the simple form at the low frequency limit / — > 1 as presented above. This normalization is just 
for illustrative purpose and not essential for our study. Note that we have correspondences such as Fj[ — > F^ and 
F% — > — F^ at order O(f ) for a plane symmetric replacement 9 — * tt/2 — 9 and <j> — * <f>. These are indeed valid at 
any order 0(f n ) (see e.g. IV. C in [HI), as we can expect from simple geometric consideration. 

With data streams A and E from a single system Bl we can make three meaningful combinations AA* , EE* and 
AE* . The expectation values for a combination C by the monopole modes ioo ancl ^oo can be written as 

<<?(/)> = [7/,c(/)Ioo(/) + 7v,c(fWoo(f)] /5. (4) 

The overlap functions 7{/ i y} i c show the relative strength of inputs ioo an d Voo to an observable (C) 0- They are 
given by the pattern functions, such as, F^(n,f) and Fj(n, /) . For example, the function ^v^E" for C = AE* is 
given by 

7V.W/) = ^J 2 dn { F a F e* - W}] • (5) 
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Here we used the definition (J5J) and equation In a same manner the function •ji^aa" is given by replacing the 
above parenthesis [• • •] with [Fj[ Fj[* + F£F£*J . The kernel [• • ■] for jv.AE'* at order 0(/ ) is given by 



I + cos 2 9 



cos9. (6) 



This factor is decomposed only with dipole (I = 1) and octupole (I = 3) patterns 0, and cannot probe the monopole 
Vbo ■ This is because responses of interferometers to incident waves have an apparent symmetry with respect to the 
detector plane, and this cancellation holds at any order 0(f n ) (see e.g. IV. C in ^3). From the same reason we 
cannot probe the mode Vbo with using self correlation, such as (AA*) or (EE*). We need independent data streams 
to capture the target Vbo- Note that the kernel for ^i^ae* becomes 8 _1 (1 — cos 2 9) 2 sin4</> at O(f ). It is written 
only with hexadecapole modes (/ = 4) 0,0], and we have ji^ae* = at 0(f°). As we see later, this is preferable 
to reduce the contamination of Too to determine the target Voo with using a combination that is a refined version of 
AE*. 

Next we consider a second system B2 in addition to the first one Bl discussed so far. With the standard configuration 
of BBO, B2 is put at position obtained b y ro tating Bl around the Z-axis by 180°, and arms of these two systems form 
a star-like shape on the XY- plane |3J, lij |lfj . But as discussed above, we can not capture Vbo with this configuration 
due to the plane symmetry of interferometers. Therefore, we study a simple case with breaking this symmetry. We 
consider to put B2 (more precisely its barycenter) on a circular orbit that has the same radius (~ 1AU) as Bl, but 
its orbital plane is inclined to that of Bl with an angle e = D/lAU < 10~ 3 in units of radian. Here the parameter D 
is the maximum distance between barycenters of Bl and B2, and its preferable scale is ~ 10 5 km, namely the same 
order as the arm-length L of BBO, as we see later. The two orbits of the barycenters intersect twice per orbital period 
^orfc(~ ly r )- in figure 2 their configurations are shown with viewing from their node. We neglect tiny misalignment 
of directions of two detector planes of order e, and only study effects caused by their relative positions. By dealing 
with the rotation of detector planes and the cartwheel motions mentioned earlier, we can follow the position of the 
B2's barycenter on the moving AYZ-coordinate attached to Bl. The trajectory of B2's barycenter (Bx, By, Bz) is 
given as 

dx = Bx/L = V3d(cosuj cos(uj + a))/2, (7) 
dy = By/L = V3d(coswsin(w + a))/2, (8) 
d z = B z /L = dcosuj/2. (9) 

Here we have defined d = D / L, and the parameter u = 2n(t/T or b) is the orbital phase of Bl around the Sun. In figure 
1 we show the trajectory of B2 for (a, d) = (0, 0.24) as dotted curves. The standard BBO configuration is recovered 
with putting d — 0. The free parameter a determines the orientation the dotted curves around the Z-axis, and we 
hereafter fix it at a = 0. 

We define two TDI modes A' and E' made from B2 system in the same way as A and E from Bl. Now two modes 
A' and E' are not on the XY-plane (except for lu = ir/2,3ir/2, ■ ■ •), and this introduces a phase shift e~ % ' znz — 
l — ifdznz + 0(f 2 ) for their pattern functions F^, E , (including information of position in the AFZ-coordinate) from 
the previous ones F^ E . When we take a combination AE'* (or almost equivalently EA'*), this phase shift generates a 
multiplier factor cos9(— n£) to eq.© at order O(f). Consequently, the combination AE'* can capture the monopole 
mode Voo at O(f), since we have a kernel proportional to (1 + cos 2 9) cos 2 9 > for circular polarization V(n). One 
the other hand, the combination AA'* (or equivalently EE'*) can be used to detect the total intensity Zoo by the 
correlation technique, as for the standard choice with d = 0. But it is important to check how the overlap function 
7j,AA'* is reduced with taking finite distances d ^ 0. 

We take a closer look at these aspects with including all the higher order effects 0(f n ). Relevant overlap functions 
are numerically evaluated, and some of the results are shown in figure 3. In this calculation we included not only 
the phase shift induced by the relative bulk positions between Bl and B2, but also the effects by the finiteness of 
the arm-length L |12(. As is well known for aligned interferometers, the function |7j aa'*| approaches 1 at the low 
frequency limit. We have the following asymptotic profiles; |tj- aa'*| = 1 + 0(f 2 ), \jv,AE'* \ = 2/|d z |/3 + 0(f 3 ) and 
1i,ae" = 0(df 2 ). The first nonvanishing term of ^i^ae 1 * is determined only by dx and dy . At d = corresponding 
to the simple traditional choice with putting Bl and B2 on the AY-plane, we have identically jv,AE'* = and the 
monopole Vbo cannot be measured, as explained earlier. When we increase the separation d (in figure 3: long-dashed 
curves — > solid curves — > short-dashed curves), the combination AA'* loses sensitivity to the total intensity loo from 
larger /. But the function \^v,ae" \ becomes larger at small / as shown by the asymptotic behavior, and the data 
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FIG. 1: Three spacecrafts (SCs) of a LISA-type system (Bl) are shown with circles that are put on the XY-plane. The 
barycenter of three SCs is at the origin (X,Y,Z) — (0,0,0). The bottom figure shows the projection of the structures to 
the XZ-plane. At low frequency limit the responses of two TDI modes A and E can be effectively regarded as those of two 
L-shaped interferometers with shown orientations. On the coordinate system (X,Y, Z) fixed to the first system Bl, the center 
of the second system B2 moves on the the dotted curves for the specific parameter choice (a, d) — (0, 0.24). A typical snapshot 
of Bl (solid line) and B2 (dashed line) is shown on the right side. The orientations of two effective L-shaped interferometers 
for A' and E' modes are also given. 




FIG. 2: Configuration of two orbital planes for Bl and B2 seen from their node (orbital phase: lu = 7r/2). 



(AE'*) get better sensitivity there to the target Vbo- In figure 3 the results for ji^ae'* are presented as a reference to 
show potential contamination of the total intensity loo for measuring the target Voo from the data (AE'*). We do not 
go into this effect. But, in many cases, it would be possible to estimate and subtract this contamination relatively 
well, as the intensity loo is generally determined better than the target Voo with using observables such as (AA'*}. 
Note also that this contamination might be somewhat reduced by adjusting free parameters including a. 



CORRELATION ANALYSIS 



We are now in a position to discuss how well we can estimate the monopole Voo of circular polarization of stochastic 
gravitational wave background. We suppose that noises of relevant data streams {A, E, A' , E'} are not correlated and 
have identical spectrum S(f) (as usually assumed for BBO). Then the signal to noise ratio (SNR) for detecting Jqo 
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FIG. 3: Overlap reduction functions for the unpolarized I and circularly polarized V modes for parameters u) = 0.5 and a = 0. 
The long-dashed curve is result with d = (traditional choice) for which we have yv,AE'* = Ii.ae 1 * = 0. The solid curves are 
for d = 1, and the short-dashed ones are for d = 5. The functions |7/,aa'*| approach 1 at / —> 0, while the functions yv,AE'* 
and Ji,ae'* have asymptotic behaviors 0(/ 1 ) and 0(f 2 ) respectively. For a — we have Ji^ae'* = (within numerical errors) 
at u> — 0. 



with using the combination AA'* is written as |]| 

^=»@) , r-'r* iTw "^r (/) ' ' 

where we have used the familiar quantity ^Igw(J) = 47r 3 / 2 / 3 /oo(/)/3i?o (Ho'- Hubble parameter) and T b s is the 
observational time. We also assumed that the amplitude of the background is much weaker than the detector noise. 
This corresponds to a situation when the correlation technique is effective. The signal to noise ratio SNRy for 
the target Vbo with using the combination AE'* is given by replacing the simple amplitude ^lGw(f) m ea. lfTU)l with 
the polarized one Htlcwif)- Here the parameter n = Vbo(/)/-^oo(/) is the polarization degree and we neglect its 
frequency dependence in this paper. 

The overlap functions 7{/,y},c(/j depend on time t through the change of the relative position between Bl and 
B2 as shown in figure 1. With the designed noise spectrum S(f) for BBO, we numerically evaluate the integrals for 
SNRi and SNRy (with n = 1) as a function of the maximum separation d = D/L between Bl and B2. We take 
the observational time T f, s as a natural number in units of the orbital period T or b, and assumed a flat spectrum 
^Gw(f) = const. In figure 4 ratios SNRi(d)/ SNRi(0) and SNRy(d)/ SNRj(0) are shown and these are our central 
results. Since we have two relevant sets {AA'* , EE'*} and {AE'* , EA'*} for measuring I 00 and Vbo respectively these 
ratios can be effectively read as the results for the total network formed by Bl and B2. When we increase the distance 
d, the sensitivity for the intensity Iqo decreases monotonically due to reduction of the overlap functions ji^aa" as 
seen in figure 3, but the sensitivity for circular polarization increases for separation d up to ~ 12. If we take d = 5, the 
ratios are SNR I (5)/SNR I (0) = 0.93 and SWi? y (5)/SWi? 7 (0) = 0.24. This means that the detection limit for the 
intensity Slew becomes slightly (~ 10%) worse compared with the simple conventional choice at d = 0, but we can get 
essentially new "sensitivity" to investigate circular polarization of gravitational wave background. For a background 
with a flat spectrum at SIgw = 10~ 15 , BBO with d — has potential to detect it at SNRi(0) = 251 by lOyr 
observation In other words its detection limit is written as tlcwMm = 2 x 10 _17 (S'A r i?//5)(r o ( );j /10yr) _1 / 2 . If 

we take d = 5 in stead of d = 0, the limit becomes H,GW,Um = 2.2 x 10 _17 (5A r i?//5)(T o f )s /10yr) _1 / 2 and we have the 
detection limit for circular polarization degree at ILu m — 0.0&(nGw/W~ 15 )~ 1 (SNR v /5)(T obs /10yr)~ 1 / 2 . 

The author would like to thank A. Taruya and T. Tanaka for comments, and A. Cooray for various supports. This 
work was funded by McCue Fund at the Center for Cosmology, UC Irvine. 
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FIG. 4: The signal to noise ratios for detecting Zoo and Voo from combinations (A A'*) and (AE'*). The results are normalized 
with SNRi(0) for the simple orbital choice d — 0. The shape of BBO noise curve is used and a flat spectrum £Igw = const is 
assumed. The polarization degree IT is set at IT = 1 for SNRv(d). 
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